Angle-dependence of the frequency correlation in random photonic media: the 
diffusive regime and its breakdown near localization. 
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The frequency correlations of light in complex photonic media are of interest as a tool for charac- 
terizing the dynamical aspects of light diffusion. We demonstrate here that the frequency correlation 
shows a pronounced angle dependence both in transmission and in reflection geometries. Using a 
broadband white light supercontinuum, this angle dependence is characterized and explained theo- 
retically by a combination of propagation effects outside the medium and coherent backscattering. 
We report a strong dependence of the coherent backscattering contribution on the scattering strength 
which cannot be explained by the diffusion theory. Our results indicate that coherent backscattering 
of the frequency correlation forms a sensitive probe for the breakdown of the diffusive regime near 
localization. 

PACS numbers: 42.25.Dd, 42. 25. Fx, 78.67.-n 



I. INTRODUCTION 

Nanostructured photonic media, such as photonic crys- 
tals, metamaterials, and random media, are of interest 
their potential of achieving nanoscale control over the 
flow and emission of light. The wave nature of light re- 
sults in a range of interference phenomena, which modify 
the propagation of light through complex photonic mate- 
rials, in analogy with the theory of electrons in solids. In- 
terference effects lead to new correlations between optical 
fields at different angles and frequencies after interaction 
with the mediumii For random or partially disordered 
photonic media, methods have appeared that allow for 
a statistical interpretation of the spectral width of the 
transmission and reflection correlations in terms of the 
distribution of dwell times of light inside the mediumi^r— 
The short-range frequency correlation, known as C^'^\ 
which is present in the spectrum after interaction with 
the medium, has evolved into a broadband tool for char- 
acterizing the dynamic transport parameters of photonic 
nanomaterials 

Here, we explore the angle-dependence of the frequency 
correlation both in transmission and reflection geome- 
tries. Previous theoretical and experimental studies of 
the C^^-* frequency correlation have not considered its 
dependence on scattering angle<^ The angle-dependence 
of the angular correlation function in transmission was 
studied by Li and Genacki, who showed that the far- 
flcld and near-field correlations form a Fourier-transform 
pair. We derive a similar relationship for the frequency 
correlation, which however is shown to produce an angle- 
dependent narrowing of the correlation width. In addi- 
tion, we investigate for the first time interference con- 
tributions to the frequency correlation in backscattering 
resulting from coherent backscattering of light. For pho- 
tonic media, time-reversal symmetry dictates that for ev- 
ery path contributing to the reflectance, a time-reversed 



path is present that interferes constructively in exact 
backscattering. This effect, which is a manifestation of 
weak localization, is known in photonics as the coherent 
backscattering (CBS) effectj^i^ The constructive interfer- 
ence in coherent backscattering results in a cone of en- 
hanced intensity versus angle around the backscattering 
direction, which is of order unity in intensity for all scat- 
tering strengths. The angular extent of the CBS cone 
increases proportional to the photonic strength (fce^s)"^ 
of the medium, where fee denotes the wavevector of light 
in the medium and (.b denotes the (Boltzmann) transport 
mean free pathJ^ 

Effects of coherent backscattering have been ob- 
served in the time-correlation function of colloidal 
suspensionsii and in the angular correlation function 
('memory effect').— However, the direct contribution 
of time-reversed light paths has never been investigated 
for the frequency correlation. A pronounced angle- 
dependence of this correlation can be expected from 
the relationship between the time duration of a short 
pulse after diffuse propagation and the spectral correla- 
tion functional, combined with the observation of angle- 
dependent changes in the pulse propagation time in 
scattering mediaJ^ We present here both a theoretical 
and experimental description of this effect, using a se- 
ries of strongly scattering scattering media with scat- 
tering strengths ranging from the diffuse scattering to 
the strongly photonic regime. By varying the scatter- 
ing strength, we study the breakdown of the diffusion 
approximation for the correlations in coherent backscat- 
tering. We report results showing that the frequency 
correlation depends on scattering strength in a new and 
previously unexplored way, and which indicate that the 
frequency correlation may provide a more sensitive means 
for accessing long light paths compared to the conven- 
tional intensity CBS effectJ^ Generally, the frequency 
correlation can be used to obtain information on the dis- 
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tribution of photon dwell times complementary to time- 
resolved experiments, which are usually not easily ex- 
tended to angle-dependent measurements.—"— In addi- 
tion, we discuss the possibility that the combination of 
CBS and frequency correlations yields additional higher- 
order correlation effects, which are unique for the four- 
field correlation and are not found in conventional inten- 
sity measurements. 



II. THEORY 



A. Angle-dependent frequency correlations in 
transmission 

Following the general definitioni, the short-range fre- 
quency correlation function C^l[u;,Q) between a single 
scattering direction u and a single incident direction a is 
given by 
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FIG. 1: (color online) Schematic picture of light transport 
through a diffusive slab, indicating injection and extraction 
planes at Zi and z"^ , Gaussian source profile (shaded area, 
red), and the additional path length difference ^ outside the 
medium. 



a function of frequency 

iJ(qx, z, z'; t)^^ I dne*f"i7(qi, z'; fi) (4) 



with tfl^a the transmission coefficient for the electric field 
and r2 the frequency difference between two fields. The 
short-range frequency correlation has been well stud- 
ied both experimentally and thcorcticallyiii^"— We follow 
here the derivation by De Boer et al. in a slightly differ- 
ent notation, where we make use of the 'intensity on a 
scatterer' I, i.e. the intensity at the last scatterer in the 
medium.— The intensity I is related to the source S 
by an intensity propagator H which contains the 4-point 
ladder vertex with attached incoming Green's functions. 
After Fourier transformation of the tranverse coordinate, 
I can be written as 



I{q±, z;t) 



dt' 



(2) 



:H{(i^,z,z';t-t')S{q^,z';t'). 



Equation [2] can be readily extended to combinations of 
fields with different frequencies by modifying the source 
S. In the following, wc consider I as a complex field- 
field correlator instead of a real-valued intensity. The 
incident light beam is described by a Gaussian intensity 
profile with waist po, which is injected into the scattering 
medium at an injection depth z^. For the frequency corre- 
lation, we consider a mixed-frequency source containing 
two oscillating fields with frequencies uj and lu + This 
yields an oscillating phase e'^* in the time response. Af- 
ter Fourier transformation, the injected source field pro- 
file, S, is given by 
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where f] (|q±p + l/DTa - D being the dif- 

fusion constant of the light. The effect of internal reflec- 
tions can be included in H following for example Zhu ct 
ali^ The diffuse propagator H in principle depends on 
the optical frequency lo through the dispersion of the pa- 
rameters VE^ £bi Ta, and D in Eq. ([S]). However, in our 
experiments we will consider these parameters constant 
as their variation with frequency in the visible spectrum 
has been shown to be smooth for the materials under 

Studji^^^>22. 

The expression for I{q±, z, t) follows from Eq. — 

as 
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where the second line is obtained by the integration over 
t' of e''^ , yielding 5{fl — fl'). Here the final expres- 
sion for I(|q^ |, z; £7) is 



S{q±,z;t) Trp^Io£BS{z - Zi)e 



(3) 



The propagator H oi a finite slab is considered, which 
will be most convenient for further use when written as 



i{q±,z;n) = 
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FIG. 2: Calculated spectral narrowing of the frequency cor- 
relation An(^)/Af2(0) against frequency correlation width 
An(0)/27r, for various values of the incident beam width po 
of 25 ^im (thick line), 250 fim (line) , and 2.5 mm (thin line). 

The transmitted intensity in the far-field can be ob- 
tained by attaching outgoing Green's functions to the 
mixed-field correlator I. In the most simple approxima- 
tion, we assume all the energy to be emitted from an 
ejection plane zj = L — 2to — with a weight of the 



ejection function given by ^b/ISttj ^^'^^ From the ejec- 
tion plane, spherical waves are emitted into the far-field. 
Since the detector collects intensity in an angular range, 
we write down the expression of the far-field intensity 
in direction u at a wavefront emitted under an angle 
with the normal. This is given by the integral of inten- 
sities over all positions rj_ on the ejection plane times 
the Green functions containing the corresponding path 
length difference ^ = ■ u = |rj^ | sin 9 cos (f>, as shown in 
Fig. [2 

i— (rn+ri -u) 

5u(ro,r_L;w) « , (8) 

47rro 

where rg denotes the distance on the optical axis from the 
detector to the center of the ejection plane (i.e. rj^ = 0) 
and is the azimuthal angle of with the radiation 
plane including 6. 

The angle-dependent far-field intensity at a distance 
To is obtained by attaching the outgoing Green's func- 
tions to the near-field intensity Eq. ([7]). This results in 
a mixed-frequency bistatic cocfficicn1(2^ in transmission 
and for perpendicular incidence, 7?'q(^): given by 



^A"^ / dr_L.gu(r-o,ri;w)gi(ro,r_L;tj + f7)i(ri,zj,rj) 

I) i — ro 



drie*<^o''^'"i(r_L,zJ;fl), 



(9) 



where A denotes the detector area following the definition 
of the bistatic coefficient)^ The second line of Eq. ^ 
can be identified as a Fourier transform of I(r_L, zj, fi) 
to a transverse momentum coordinate ~"^'"^ u i , where 

Co ^ ' 

sin 9 is the projection of the scatting vector u onto 
the ejection plane. This Fourier relationship between the 
near-field and the far field spectral correlation follows the 
van Cittcrt-Zernike theorem and is an example of a more 
general theorem dealing with spectral invariance of light 
on propagation^. A similar Fourier relationship holds 
for the angular correlation function, as was shown by Li 
and Genackii. The Fourier transform results in a final 
expression for the scattering bistatic coefficient 
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with X the near-field correlator in momentum space, 
Eq. (O. The correlation function of Eq. ([!]) can be 

calculated by identifying that o('^ + a^'^-') ^ 

^fl o^^)' which the correlation in transmission fol- 



(11) 



Equation ([T0|) shows that the far-field frequency cor- 
relation gains an additional angle-dependent dcphasing 
which can be written as a transverse momentum which 
is proportional to the frequency difference fi/co and the 
angle 9. 

The spectral correlation of Eq. [TT] has a characteris- 
tic spectral width Ar2(0), which in the forward direction 
is approximately given by Af2(0)/27r ~ lAQD/Li^J^ For 
nonzero angles, the spectral width is reduced, as is shown 
in Fig. [5] where we have calculated the spectral narrow- 
ing Af2(7r/2)/Afi(0) as a function of the full-width at 
half maximum of the spectral correlation function. Fig- 
ure[2]shows that the spectral narrowing depends strongly 
on the combination of the incident beam size po and the 
initial spectral width. For typical beam sizes of around 
25 — 250 /xm, the narrowing effect is pronounced for cor- 
relation widths in the range 10 — 100 cm~^. This fre- 
quency width corresponds to dwell times Td — Tr/AfJ 
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FIG. 3: (color online) (a) Four-field diagrams contributing 
to the frequency correlation C^i^\lo,Q,), from left to right 
LL*, CC*, LC*, and CL*. (b) Frequency full-width at half- 
maximum AQ, normalized to width at backscattering ^ = 0, 
for the different diagram contributions, (c) Relative contri- 
bution of various diagrams to the total C'^' correlation, (d) 
Total calculated broadening of correlation function against 
backscattering angle. 



in the range 0.1-1 ps, or several hundreds to thousands 
of scattering events. This regime applies to many types 
of photonic nanomaterials with thickness in the range of 
hundreds of nanometers to several micrometers, demon- 
strating that the angle dependence of the frequency cor- 
relation in Eq. PT|) is a sizeable effect in many situations. 



B. Angle-dependent frequency correlations in 
reflection 

While the C^^-* frequency correlations in transmission 
are governed by long light paths scaling with the slab 
thickness L, the correlation in reflection is governed by 
paths of the order of the mean free path Ib- In addition, 
the reflection geometry is characterized by the presence 
of weak localization in the backscattering direction. 

We start again from the frequency correlation of 
Eq. ([T]), where this time we consider reflection coeffi- 
cients Q for the electric field. For the case of reflec- 
tion, we expand the numerator of Eq. ([1]), analogous 
to Berkovits and Kavebi^, into combinations of diffu- 
sive and maximally-crossed multiple scattering contri- 
butions, respectively denoted here by L (diffuson) and 
C (Cooperon) [cf. Fig. [Ha)].— The spectral correla- 
tion function consists of four diagrams, consisting of two 
diffusons (LL*), two Cooperons (CC*), and two mixed- 
diagram contributions (LC* and CL*). We write the 



correlation as the sum of two mixed-frequency bistatic 
coefScients 7!. -{uj, V.) and 7^ -{uj, V.) 



(12) 



where is a normalization factor given by = 
37fi a(^' 0)"^. Analogous to the derivation in transmis- 
sion, the bistatic coefficient for the ladder contribution 
can be written as 



sin 6 

Co 



Ui,zf;fl), (13) 



where = |£ denotes the ejection plane in reflection, 
based on the average and the weight of the ejection func- 
tion as defined in Ref. [l^. For the case of the coherent 
backscattering contribution, the solution is obtained by 
adding the incident and outgoing plane waves on posi- 
tions r I , r'l — 



7c(q±,w,l7) 



d^r'i d^r le-'i^-f'-^-'-'^) 



x7?(ri -rl,z„ze,f^)^(ri,z,;l)), (14) 



where = fcou± sin 9 denotes the transverse wavevector 
difference of the incident and outgoing waves. Following 
the derivation of Sec. Ill Al we find a similar expression for 
the bistatic coefficient with the addition of a wavevector 
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The angle-dependent frequency width of the various 
diagrammatic contribution to the frequency correlation 
is shown in Fig.|3Ib), while their relative amplitude con- 
tributions to the total correlation is shown in Fig. |3Uc). 
In exact backscattering [9 — 0), 7c equals 7;, and the con- 
tributions of the four combinations LL*, CC*, CL*, and 
LC* are equal. For nonzero wavevector qx , the CBS con- 
tributions broaden spectrally due to the reduced weight 
of long light paths [cf. Fig. EJb)]. However, the CBS 
contribution decreases rapidly in intensity for large scat- 
tering angles [cf. Fig. IHtc)], resulting in recovery of the 
purely diffusive correlation at large angles. The resulting 
frequency broadening of the total correlation, including 
the contributions of all diagrams, is shown in Fig. [3l^d). 

The angle dependence in the coherent backscatter- 
ing frequency correlation shows similarities to the time- 
correlation function of weakly scattering colloidal sus- 
pensions in backscattering] ^^'^^ These effects are related 
because of the presence of the same path length distribu- 
tion in the dephasing caused by moving particles. How- 
ever, whereas the time-correlation is governed by the slow 
movement of scattering colloids, the frequency correla- 
tion in our case directly relates to the photon transport 
time in the random medium with quenched disorder. 
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FIG. 4: (color online) Experimental setup for white-light, 
angle-dependent frequency correlations in transmission and 
backscattering geometries. Dashed lines indicate components 
which are only present in backscattering configuration. 



III. METHODS 

Correlations were measured in the frequency domain 
using a broadband technique developed by us as de- 
scribed in Ref0. Spectra containing individual frequency 
speckles were collected as a function of scattering angle 
using a rotation stage and a fibcr-couplcd grating spec- 
trometer, as shown in Fig. 2] Polarization filters were 
used to select linear polarization channels for the inci- 
dent and detected light. In the coherent backscattering 
experiments, a quarter wavelength plate was placed di- 
rectly in front of the sample to select circular polariza- 
tion channels. Circular polarization is used to suppress 
single-scattering contributions in the coherent backscat- 
tering cone, which would otherwise result in a reduced 
enhancement factorial The spectral correlation function 
was obtained at each angle by averaging the correlation 
functions from 50 spectra taken at different sample po- 
sitions. The correlation function was averaged over a 
spectral range from 550-850 nm, yielding an average over 
several thousands of individual speckles. The configura- 
tion could be switched between a transmission geometry 
and a backscattering geometry as illustrated in Fig. |4l 
The backscattering setup had an angular resolution of 
1.2 mrad. For large-angle backscattering measurements, 
the beam-splitter was replaced by a small reflecting prism 
following Ref. I2TI . For the transmission experiments, the 
beam waist at the sample position was varied by using a 
set of lenses of focal lengths ranging between 5 cm and 
60 cm. The resulting beam waist at the sample posi- 
tion was determined using a knife edge method, yielding 
values of po between 26.4 ± 0.7 /im and 328 ± 16 pm. 

For the transmission measurements in Sec. IIV Al we 
used two thin slabs of Ti02 powder with thicknesses of 
3.4 fini and 6.3 /im. The mean free path of the sam- 
ple ranges between 0.55-0.9 ^m in the visible part of the 
spectrumi^ The thin Ti02 slabs have been investigated 
before using time-resolved transmission, and were found 
to be within the limits of classical diffusion description.'^S 
In the reflection experiments of Sec. IIVBI we study a 
series of random scattering media with different inverse 
photonic strengths keis, covering nearly an order of mag- 
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FIG. 5: (color online) (a) Experimental frequency correlations 
after transmission through a 3.4-^m slab of Ti02 powder, for 
a few typical forward angles, (b) Experimental correlation 
width AQ/2n against transmission angle for a beam waist 
Po of 26.4 ± 0.7 p.m (dots, black), 55 ± 2 pm (circles, red), 
120 ± 4 nm (upright triangles, green), 328 ± 16 fim (open 
triangles, blue). The three black open circles with arrows 
correspond to the experiments displayed in (a). 



nitude from fc^^B = 3.6 ± 0.8 to kjB = 26 ± 3. The 
materials under study include samples of etched porous 
GaP which have been under investigation earlier for the 
possible effects of localizationji^iii These samples have 
been shown to be of sufficient thickness and with a low 
residual absorption for sensitive measurements of long 
light paths. In addition, we also include in our stud- 
ies two layers of GaP nanowires, which have been shown 
to be among the most strongly scattering nanomaterials 
available today<2^ 



IV. EXPERIMENTS 

A. Frequency correlations in transmission 

We investigated experimentally the angle-dependence 
of the frequency correlation in transmission for various 
values of the beam waist po- Figure [5ja) shows typi- 
cal spectral correlation functions for po = 328 ± 16 pm 
at transmission angles of 0.04, 0.7, and 1.4 rad, 0.0 rad 
indicating the forward scattering direction. The exact 
forward direction was not measured as this contained 
contributions from the coherent beam for the optically 
thin samples {L/£b — 5 for the 3.4 pm thick slab). The 
correlation functions were normalized to the second data 
point, to remove any uncorrelated noise which is accu- 
mulated at zero frequency shift. A pronounced narrow- 
ing of the spectral correlation is observed for increasing 
transmission angles. This narrowing is further character- 
ized in Fig. [Sjb) where the fuU-width-at-half-maximum 
An of the spectral correlation is plotted against trans- 
mission angle, and in Fig. [5] showing the narrowing at 
1.4 rad for the two different slabs. Clearly, the narrow- 
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FIG. 6: (color online) (a) Experimental values of the spectral 
correlation width ASl at 1.4 rad, normalized to the width at 
rad, against beam waist po, for Ti02 slabs of 3.4 /im (cir- 
cles, red) and 6.3 fim (triangles, blue) thickness, (lines) Cal- 
culated spectral narrowing of the Cy^' correlation of Eq. (|lip 
for transmission. 



ing depends on beam waist po, in agreement with our 
theoretical calculations using Eq. pT|) . as illustrated by 
the hues in Fig. El For the calculations we used values of 
L, £, and D obtained from other experiments and which 
produced the correct width of the frequency correlation 
function at zero angle. Therefore, the curves in Fig. [6] 
do not contain any adjustable parameters. The increase 
of the slab thickness from 3.4 ^m to 6.3 ^m results in a 
longer dwell time and thus in a reduction of the correla- 
tion width. This results in a correspondingly less strong 
angle-dependence, following Fig. [21 We point out that 
the agreement with theory is limited to beam sizes that 
are smaller than the diameter of our detector, which in 
our experiments is set by the 500 fim diameter fiber. We 
conclude that the model including angle-dependent prop- 
agation outside the medium, as given by Eq. (jlip . pro- 
vides a good quantitative description of our experimental 
results, including the scaling of the effect with po- 



B. Frequency correlations in coherent 
backscattering 

In the following, we examine the single-channel fre- 
quency correlation around the coherent backscattering 
direction. Compared to the transmission experiments, 
the reflection measurements require slabs of sufficient 
thickness and scattering strength to resolve the coher- 
ent backscattering effect. We have chosen to study a 
series of slabs of porous GaP with different scattering 
strengths, as these have been shown to cover the tran- 
sition from diffuse scattering to the strong scattering 
regime near localization . ^^'^^ In addition, two strongly 
scattering porous GaP slabs were measured before and 
after infiltration with 1-dodecanol. Infiltration with 1- 
dodecanol lowers the photonic strength by a factor three; 




Frequency Ql2n (cm"^) Frequency Ql2n (cm^) 



FIG. 7: (color online) Experimental correlation functions in 
reflection (normalized) for porous GaP samples before (a) and 
after (b) infiltration with 1-dodecanol. Symbols indicate cor- 
relation in backscattering (diamonds, black ai 6 = rad) 
and at the maximum broadening (dots, red at respectively 
9 — —0.06 rad and 6 = 0.024 rad). (ii), (iii) correspond to 
Fig.H 

this method allows unambiguous separation of the ef- 
fects of scattering strength from other, sample depen- 
dent contributions!^ Figure [71 shows the experimental 
correlation functions of a porous GaP sample before (a) 
and after (b) infiltration, and for angles corresponding 
to exact backscattering {6 = rad) and angles where the 
largest change in spectral width occurs. The frequency 
correlation functions in reflection show a triangular peak 
which is slightly rounded by our spectral resolution of 
around 3 cm~^. Upon infiltration with 1-dodecanol, the 
frequency correlation narrows indicating an increase of 
the photon dwell times resulting from the increase in 
mean free path. 

In the same experimental configuration as used for the 
frequency correlations, conventional intensity CBS cones 
were obtained by collecting the average intensity. Fig- 
ure |5Ja,b) shows a typical combination of intensity CBS 
cone (a) and the 1 /e frequency half- width of the spectral 
correlation function (b) as a function of backscattering 
angle. The black dots represent measurements taken in 
the helicity conserving channel, while the open triangles 
(red) in Fig. jSl^b) are measurements in the polarization 
nonconserving channel. The narrowing of the spectral 
correlation at large angles can be attributed to the de- 
phasing contribution of Eq. ([TO)) . As the path length dis- 
tribution in reflection is governed by short light paths of 
the order of the mean free path £b, this angle-dependent 
narrowing stays prominent even when L increases to in- 
finity, i.e. for semi-infinite slabs. 

An additional feature is observed in the helicity con- 
serving channel which resembles the calculated response 
of Fig. |2Id). For the narrow angular range around 
backscattering where the effects of the CBS cone are 
most pronounced, the additional narrowing due to de- 
phasing is relatively small. We therefore analyze the 
broadening associated with the coherent backscattering 
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FIG. 8; (color online) (a) Experimental intensity CBS-cone 
measured for a porous GaP slab, with (line) fit for keis = 
4.7 ± 0.6. (b) Frequency correlation 1/e half-width Af2/27r 
against backscattering angle for circularly polarized conserv- 
ing (dots, black) and linearly polarized nonconserving (red, 
open triangles) channels. 



effect by defining an experimental spectral broadening 
parameter Af2exp(^) as the ratio of the helicity conserv- 
ing and nonconserving widths AJ7 at every angle 9. Re- 
sults are shown in Fig. [9la,b) for 4 different samples 
with scattering strength ranging from keis = 26 ± 3 to 
keis = 3.6 ± 0.8. The wavevector k^ includes the effec- 
tive refractive index of the scattering medium, ng, rang- 
ing between 1.4 and 2.0 for the materials under study. 
Backscattering corresponds to 6* = 0.0 rad. Conventional 
CBS cone (i) corresponds to the most strongly scattering 
nanowirc material of Ref. [2^ . CBS cones (ii) and (iii) in 
Fig. HJa) were taken from the same porous GaP sample 
before and after infiltration with 1-dodecanol, consistent 
with earlier measurements on the same samplei^^ Cone 
(iv) corresponds to a less strongly scattering porous GaP 
layer. The red lines in Fig. IHJa) show fits using a finite- 
slab model^ including an internal reflection correction, 
yielding values of keis as indicated in the figure. All 
conventional CBS cone shapes can be described by the 
finite sample thickness or limited angular resolution, the 
estimated absorption lengths exceed the sample sizej^^ 

The characteristic anglc-dcpendcncc of Aficxp is found 
in Fig. [DJb) for all samples under study, however with 
different magnitudes of the broadening effect. The 
angle-dependent broadening can be reproduced theoreti- 
cally using Eq. (|12p using the known sample parameters 
(dashed lines, blue). This model gives good quantitative 
agreement for the least strongly scattering material with 
keis = 26. For the more strongly scattering samples with 
kgis < 10, the magnitude of the experimental broaden- 
ing is significantly smaller than would be expected on the 
basis of the C'-^' theory. In our model analysis, we take 
into account the effect of known finite slab thicknesses, 
which cuts off of light paths above a critical length in the 
propagator of Eq. ([5]). This correction is especially im- 
portant for the nanowirc material (i) for which L/£ ~ 25. 

The red lines in Fig.|njb) correspond to the C^^-* theory 
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FIG. 9: (color online) Left (a) Conventional intensity-CBS 
cones for different random scattering media, with from top 
to bottom: compressed nanowires, photoanodically etched 
porous GaP in air and in dodecanol, and etched porous GaP, 
with values keis ranging from 3.6±0.8 to 26±3. (Lines, red) 
Theoretical fits of conventional CBS^. Right (b) Experimen- 
tally determined increase in the width of the spectral correla- 
tion function Aficxp. (Lines) Theoretical calculations of using 
Eq. (|12p for known material parameters (Aficaic, dashed lines, 
blue) , and including deviation parameter S characterizing the 
long range correlations (Afioxp, lines, red). Vertical bar de- 
notes deviation from theoretical maximum [S = 0, horizontal 
line, blue). 



with an amplitude fitted to match the experimental data. 
We find that the theory can be scaled to our experimen- 
tal data to good agreement using only a single, angle- 
independent scale factor. Therefore, in order to quantify 
disagreement with the diffusion model, we introduce an 
angle-independent parameter S describing the percent- 
age of deviation from diffusion, which is defined as the 
ratio of experimental to calculated broadening curves, 
offset by the baseline of unity, i.e. 



S = 



(16) 



The finite values of the slab thickness results in reduced 
amplitudes of the calculated curves, as can be observed in 
Fig- IH^) -The vertical bars (blue) in Fig.[51Jb) indicate the 
deviation of the experimental results from the C'-^' the- 
ory corrected for slab thickness (S ~ 0, blue lines). The 
resulting deviation from classical diffusion S is presented 
for all studied samples in Fig. [TUTa) as a function of the 
inverse photonic strength keis- The experimental data 
show a strong increase of S toward fce^s ~ 3.6, which is 
unaccounted for by finite sample thickness and absorp- 
tion, as can be concluded from the CBS-intensity cones 
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FIG. 10: (a) Experimentally observed deviation from the sim- 
ple C'^^ CBS model, S, against keis- (Open symbols) Pairs 
of scattering samples before and after infiltration with do- 
decanol. Lines represent C'^-* value of S = (dash-dotted 
line) and Eq. S = ^{kdB)~'^ (solid line), (b) Spectral 

1/e half width Afi of the C'^' correlation function at exact 
backscattering 6 = Q against kels for the different materi- 
als under study. Line denotes fit using the classical diffusion 
theory with Af7 ~ D/e.% « 1/fcefs. 



of Fig. m^a). Infiltration with dodecanol significantly re- 
duces the deviation from classical diffusion, as indicated 
by the open symbols in Fig. [TO] Each type (hollow cir- 
cles, diamonds) represents one particular sample before 
and after infiltration. The pronounced effect of infiltra- 
tion agrees with our expectation for an effect depending 
on the scattering strength, rather than on other sample- 
dependent properties. The frequency width AJ7 of the 
correlation in exact backscattering, shown in Fig. [TUT b). 
follows a l/keig dependence which can be interpreted 
from classical diffusion in the following way. For clas- 
sical diffusion, the frequency width Afl is proportional 
to D/£^, which in combination with D ~ ve£bI'A yields 
the observed inverse linear dependence on Ib- A renor- 
malization of the diffusion constant^^^ would lead to a 
narrowing of the frequency correlation, which is not ob- 
served in our experiment. 



V. DISCUSSION 

We interpret the general trend in Fig. [TOla) starting 
from the diffuse scattering regime. In the CBS cone, the 
wings of the cone correspond to very short transverse 
distances between the point of entry and the point of 
exit of light in the medium, of the order of the mean 
free path Ib'^ The frequency broadening Afioxp in the 
wings is sensitive to the phase delay of these short paths 
relative to that of the total distribution. In the C'^-' ap- 
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FIG. 11: Higher-order diagrams in the frequency correlation 
in reflection, including time-reversal symmetry^S. L and C 
are ladder and Cooperon vertices, H indicates the Hikami- 
box vertex consisting of three terms.— 



proximation, the maximum frequency broadening AOcaic 
is independent of k^tB, as indicated by the dash-dotted 
line in Fig. llOf a). The deviation of our experimen- 
tal data from this line can be caused by two possible 
types of effects, which cither affect the dwell time, or 
which add new contributions to the spectral correla- 
tion. To the first type belongs the cutoff of long light 
paths and the associated distribution of decay times near 
localization! ^"'^'•^'^ The second type includes higher order 
contributions to the correlation function known as C*^^^ X 
Several authors have considered this contribution in a re- 
flection geometry] i2,27- _30 r^j^^ angular C^"^^ correlation 
is inversely proportional to the illumination area, and 
is expected to be negligibly small for the geometry of 
our experiment. In contrast, the near-field correlation in 
refiection is a local effect on a length scale of order ^, 
which does not vanish for plane wave illumination. This 
long-range near-field correlation has been calculated by 
Stephen and Czwilich as21 



27 



(17) 



Equation (flT)) shows that, for short the distances R ~ £b, 
contributing to the wings of the CBS cone, the near- 
field intensities are much stronger correlated than for dis- 
tances R» £b contributing to the center of the cone. 

At this stage, a full theoretical framework of the C^^^ 
contribution to the coherent backscattering frequency 
correlation has yet to be developed. Figure [TT] shows 
the higher-order diagrams that have to be included in 
such an extended model, following Ref. [2^. The four- 
field correlation takes place through the Hikami-box ver- 
tex, which is not present in the conventional intensity 
CBS. In addition, interference with an unscattered, re- 
duced incident intensity may be important.— In a dif- 
ferent context, these 'shadow' terms have been identified 
for their role in energy conservation in the CBS-conci^ 
Althogether, the frequency correlation CBS may contain 
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sizeable corrections to diffuse transport for values of Mb 
larger than 1 (respectively kis — 2.5 and Mb — 1-8 
in Refs. [l^ and IM ). Self-consistent theories of time- 
dependent diffusionHi^ predict a reduction of the dif- 
fusion constant in this regime, which should lead to a 
narrowing of the spectral correlation. This is not ob- 
served experimentally in Fig. Iior b). However, we point 
out that both the horizontal axis depending on £b and 
the vertical axis depending on D should be renormalized 
in this regime, making it hard to observe these renormal- 
izations experimentally as there is no absolute reference 
against which these can be calibrated. 



FIG. 12: Comparison between the sensitivity of the parameter 
S and the CBS cone roundingi^ AOb. to the suppression of 
long light paths, represented by the cutoff length L^. 



new information compared to the conventional intensity 
CBS, where C^^^ is only a backgroundi^ This is impor- 
tant as our measurements show a collapse of the CBS ef- 
fect in the correlation function, while the intensity CBS 
cone is still intact. A crucial point is that self- interference 
due to the intersection of paths is associated with a fi- 
nite intersection volume of order X^Ib, with A the optical 
wavelengthiiS This is especially relevant in the wings of 
the CBS cone, where the transverse distance between the 
path ends is only one mean free path £b- Following the 
scaling of Eq. [T71 the condition C^^^(^b) = 1 is reached 
when kjB (27/2)^/2 ^ 3^7^ ^he breakdown of the 
weak-scattering limit for the C*-^^ correlation is thus ex- 
pected to occur far above the localization limit ke£B ^ 1 
in the wings of the CBS cone. The line in Fig. [TOla) il- 
lustrates this scaling, where we hypothesized a deviation 
S of zero in the diffuse limit up to full breakdown S = 1 
of the C*-^-* model due to path intersection. 

Apart from possible higher-order contributions in the 
spectral correlation, wc consider the sensitivity of the 
CBS-induced broadening of the frequency correlation to 
changes in the path length distribution. In particular, we 
compare the new technique to measurements of the CBS- 
cone rounding in intensity measurements Figure [12] 
shows both the cone rounding angle A^^^ and the sup- 
pression of the frequency broadening S, as a function 
of the path length cutoff defined by a (total) length L^. 
For simplicity we choose an exponential cutoff such as is 
the case for absorption in Eq. ([5]), however alternative 
choices of the path length scaling are available to de- 
scribe other effects such as localization]^ Experiments on 
the same porous GaP samples have identified an anoma- 
lous cone rounding A9 of around 4 mradJ^ In Fig. [T^] 
this value corresponds to a value L^/£ of 3.7 x 10'^, or 
a cutoff length of around 1 mm for the porous GaP 
sample. In comparison, for the same L^/£ the param- 
eter S gives a value of 0.28. Experimentally we have 
observed, for the same sample as used in RefJ^, an S- 
value of 0.65 ± 0.04. Recent experiments have reported 



VI. CONCLUSIONS 

In conclusion, wc have studied both theoretically and 
experimentally the angle-dependence of frequency corre- 
lations in random photonic media both in transmission 
and in reflection. A narrowing of the spectral correlation 
was found which depends on both the scattering angle 
and the incident beam waist, which could be identified 
as an additional dephasing of the multiple scattered light 
outside the medium. This contribution has not been re- 
ported in other work on frequency correlations fir— We 
note that earlier experimental studies dealt with much 
thicker slabs for which the width of the frequency cor- 
relation was too small to observe the effect according 
to Fig. [2] However, the angle-dependent narrowing is 
sizeable for frequency widths above 10 cm~^, which cor- 
responds to dwell times below 1 ps. Thus the effect is 
relevant for many nanophotonic materials and thin films 
in which light trapping plays a role, and particularly in 
reflection geometries where the path length distibution 
is governed by short paths on the scale of the mean free 
path £b, even for very thick samples. 

In backscattering, we have found a contribution to 
the frequency correlation which can be attributed to 
the angle-dependent path length distribution in coher- 
ent backscattering. Our diffusion model is confirmed by 
experiments for scattering materials with inverse pho- 
tonic strength fce£_B — 26, however the magnitude of 
the broadening effect appears reduced for samples with 
larger scattering strength. We emphasize that a sim- 
ilar strong dependence is not found for the intensity 
CBS cones obtained in the same experiment. Thus, 
the frequency correlation in backscattering may form 
a new instrument for accessing the breakdown of dif- 
fusion theory in the limit of very strong multiple scat- 
tering. Earlier work on similar (and some of the same) 
strongly scattering materials revealed no or only weak 
corrections to classical diffusiouj^i It has been pointed 
out that time-resolved transmission is not expected to 
yield large modifications for the typical samples under 
stud}*^. The angle-dependent frequency correlation pro- 
vides a means for accessing the path length distribu- 
tion in reflection on much shorter time scales than can 
be accessed using ultrafast laser pulses. The frequency- 
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resolved technique thus yields information complimen- 
tary to picosecond time-resolved experiments which are 
usually more involved and not easily extended to study- 
ing angle-dependent phenomenaJ^ii^"— The frequency 
correlations in coherent backscattering have the potential 
to be an even more precise probe of changes in the path 
length distribution than cone rounding measurements us- 
ing the intensity CBS, making it an important tool for 
studying the strong scattering regime near localization. 
Our results call for new theoretical models for frequency 
correlations in combination with coherent backscattering 



in the strong scattering limit. 
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